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We report our exploratory study on the matching condition of HQET with QCD including 1 /m 
corrections. We introduce a new observable from the dependence of the heavy-light effective 
energy on the twisted boundary condition parameter 0, which could be used to match the kinetic 
term ly /2m. Carrying out quenched QCD simulations for fixed lattice spacing in small volumes 
with 0(a)-improved Wilson fermions, we study the 1/m dependence of this observable, from 
which the static limit and 1 / m coefficient can be extracted. We also compare our preliminary 
result with HQET. 
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1. Introduction 

One of the primary goals in particle physics is to determine fundamental parameters of the 
standard model such as the quark masses and the Cabibbo-Kobayashi-Maskawa matrix elements 
(e.g. Vtd; V'ts, Vtb) in order to test the standard model and find a clue to the physics beyond the 
standard model. However for this determination in addition to the experimental inputs from CLEO- 
c BaBar [|] and BELLE we need non-perturbative evaluation of hadronic observables 
from the first principles of QCD. Above all, simulation of lattice QCD is one of the most effective 
approaches. 

In lattice QCD, the typical lattice cutoff accessible in the present computers is t?^^ f» 1 ~ 
3GeV. But in the case of heavy quarks such as the bottom quark whose mass is m 4GeV the 
quark mass in lattice unit is so large ma 1 ~ 4 that the conventional lattice fermion action has no 
control over the discretization errors. In order to avoid this problem which has been raised of the 
heavy-quark calculation, we must take the following approaches: 

(a) conventional fermions with much finer lattice spacing, 

(b) non-relativistic QCD and their variations, 

(c) anisotropic lattices [Q, §], 

(d) heavy-quark effective theory (HQET) |§ 

However (||) is still not practical in high precision computation due to computational cost and (^ 
has theoretical uncertainty owing to the perturbative matching, while (||) and (^ have potential of 
providing formalism for high precision computation. In this report, we will give our exploratory 
study on (|) HQET. 

HQET is a theory that reproduces the low energy mode of a heavy quark in the static limit 



|@, |10p, where the action is a systematic expansion in inverse powers of the mass 

Shqet = a^I l^statW + I ^^"K^) I , (1-1) 
^statW = %{V*o + 5m)Yi.{x), =I«?'^^-^'k^), (1-2) 

5«.=^fiy«j')=i+^fi,y «(^)=i+^fi,y (1.3) 

\m J m \m'^ J m \m'^ J 

^['\x) = % Yh, = % (-^D') V^h. (1.4) 

If one treats in 1/m corrections perturbatively to a fixed order, the theory is a renormalizable theory 
with well-defined continuum limit: We presume it can accept up until next-to-leading order terms 
of 1/m and can have 1 ~ 2% accuracy. However HQET is an effective theory with unknown 
coefficients for 1/m connections terms and the ultraviolet behavior is essentially different from 
QCD. Therefore, one has to match HQET with QCD by matching conditions defined by a set of 
physical observables 

(^qcd(1/z)) = (^hqet(1/z)) , z = ML » 1 (^;an observable), (1.5) 
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where M is a renormalization group invariant mass [[T^]. This can practically be done in a small 
volume L « 0.2fm where we take the small lattice spacing within reasonable numerical cost: 
rn\,a <C 1. Then by computing the step-scaling function [[l0|], we can evolve the lattice HQET 
into coarser lattice where physical observables can be computed in large volume. In fact these two 
steps have been caiTied out for Sm^ and the axial current A^*^^ = A^^^^ in the static limit, from 
which the following results are obtained 



Fj^';^' = 253±45MeV m, mb = 4.12 ±0.07 ±0.04MeV [|T0|. (1.6) 



For a high precision calculation, non-perturbative calculation of the 1/m correction is required. 
For this purpose it is necessary to match the coefficients of the 1/m connections in the action or 
operators to begin with. 

The matching of HQET with QCD is composed of three steps: 
(e) to execute the matching condition ( |1.5| ), 



(f) to match coefficients of 1 /m-correction terms from Eq. (1.5), 

(g) to evaluate the step scaling functions. 

In this report we give study on (|) as the first step for (|) and also computing the static limit; 
especially on the search for new observables which are efficient for the determination of 1/m term. 

2. Observables 



In Ref. [13], Heitger et al. propose using two-point conelation functions in small finite volume 



as the matching condition (|1.5|). Thus we can evaluate some HQET pai^ameters by combining 
observables in varied kinematical conditions and execute heavy quark lattice simulations of both 
QCD and HQET in Lq « 0.2fm. 

Two-point correlation functions are defined by 

/a(xo) = -\ j ^/V^/'z(AoCb(y)75Ci(z)> , ^v(xo) = -\L j dV^{vMy)yM^)) , (2.1) 

where A^(x) = V/i(x)7;i75 Vb W, V'^W = %{x)ynWb{x), and C,, C, ai-e the boundary fields atxo = 0. 
Defining effective energies as 

Ef^{xo, d) = -^ln[/A(xo)], Ef{xo,d) = -^ln[^v(^o)], (2.2) 
we can determine various observables in kinematical conditions such as 



r„ = -i^E}^,y-j+3Et,-i^-Jj.AT = El-y-j-E}^,y-), (2.3, 

1 ' -eff M \ , o E-eff If _„ff /T\ _ _„ff / T 



to assess 5m, (ol^\ and Wj^"* respectively, defined in Eq. (|L^. 

It is known that the kinetic term of 1/m correction (1)2^ is difficult to evaluate it with good 
accuracy, because of the cancellation between power divergences of effective energies at r/2 and 
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r/4, which causes large systematic eiTors. Therefore it is fruitful to look for alternative observables 
for consistency checks and pursuit of more efficiency. In this work we propose an alternative 
observable: 

S-- = r,,(l/z)|e=i.o-rav(l/z)|0=o.5, (2.5) 
where 6 is defined by twisted boundary conditions 

VA(x + k) =e'VW, vix + kL) = Yix)e-"^ (it =1,2,3, T=L, 6=0.5,1.0). (2.6) 

3. Simulation methods 

Simulations are carried out both for QCD and HQET on a quenched 12"* lattice with the stan- 
dai^d plaquette gauge action at j3 = 6/gQ = 7.4802. We take the Schrodinger functional bound- 
aries [0, [T5|] with C = C = 0. The SF-boundary coefficients are c^"'°°p = 1 - 0.089gg - 0.030g^ 



[16]. 256 configurations are accumulated by the pseudo-heat-bath algorithm each separated by 200 
Monte Cai^lo sweeps. Errors are estimated by the standard jackknife method. 

In QCD, we use O (a) -improved Wilson fermion with the nonperturbative clover coefficient 
csw [ p^ and the SF boundary conditions P^y{^)\xo=o = p, P-Y{^)\xo=t = [|T8|]. Six values 
of hopping parameters are taken for the heavy quark which correspond to z = ML = 3.0, 3.8, 
5.15, 6.0, 6.6, 9.0 whereas the hopping parameter for the light quark is set to the critical value 
(k^ = 0.133961) following Ref. 

In HQET, we used the static-limit action ^stat {x) ( |1.2| ) with the boundaries for the heavy quark. 
In order to improve the numerical precision, we adopted gauge fields in Vq of ^stat {x) are smeared 



by the HYP links ||20|] with (tti, a2, "3) = (1.0, 1.0, 0.5) |gl|, The tight quark is the same as 
in QCD. 

4. Results 

Fig. |I| shows the average effective energies Pay (2.3) at 6 = 0.5, 1.0 with matching condition 



LFav C^^ss (M/Ams )z + ^( ( 1 /z)° ) in , (4.1) 



where Cmass is the matching function between QCD and HQET [ |13| ] and M/A-^ is given in 
Ref. ^n\\. We reproduced consistent result with 6 = 0.5 at j3 = 7.4802 in Ref. [p^]. As is obvious 
from Fig. |, the data for = 0.5 at j8 = 7.4802 deviate only by 2 ~ 3% for z = 3.0, 3.8, 5.15, 
whereas for z = ML = 9.0 the deviation is of 10% order due to the discretization error. Although 
there is no scaling study, we find the 1 /z dependence of the result for = 1 .0 also shows quali- 
tatively similar behavior. In fact, our data at z = 3.0, 3.8, 5.15 seem to approach the static limit 



value predicted as in Eq. (4.1 ). On the other hand, the data for z = 9.0 seem to deviate from the 
1/m scaling significantly. There is need of a stringent test in further studies, since this deviation is 
probably due to discretization enors. 

Fig. shows the l/z dependence of the observable S"™. Owing to the reparameterization 



invariance |23], the renormalization factor is not needed for this observable. Assuming that the 
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Figure 1: The averages of the effective energy at 
9 = 0.5 (drcle symbol) and = 1.0 (diamond). 
The cross symbols and the dashed line denote the 



data for = 0.5 at j3 = 7.4802 [[T9| and in the 
continuum limit [Oj respectively. 
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Figure 2: 1 jz dependence of the new observable 
^2 new rpj^g circle symbols denote our data for 
QCD and the diamond denotes the HQET result. 
The dashed line is the linear fit using three lighter 
points. 



data for z = 3.0, 3.8, 5.15 are very close to the continuum limit, we made a linear fit in 1/z 
as LE"''* = ^0 + ^1 jz- Our preliminary results are ao = 0.750(12) and a\ = 0.479(47), wliere the 
errors are statistical only. The data for z = 9.0 again deviate from the 1 /z scaling, which is probably 
due to the discretization error. The extrapolation of QCD results does approach to the HQET 
result towards the static limit, although it gives higher values than the HQET result at j3 = 7.4802. 
Further study is needed to see if this discrepancy vanishes in the continuum limit. Nevertheless, it 
is encouraging that this observable has clear 1/z dependence so that it has sensitivity to determine 
the coefficient of 1 /m-correction terms COj^' in HQET. 



5. Conclusions 

We have proposed the new observable E"*^* defined as the difference of the effective energy 
at 6 = 0.5, 1.0 for extracting the coefficients of the kinetic term in HQET. Our simulation at 
j3 = 7.4802 suggests that it has qualitatively correct HQET scahng and also good sensitivity to 
1 /m-correction term. 

Acknowledgements: The simulation has been done at the Yukawa Institute Computer Facil- 
ity and on a supercomputer (NEC SX-5) at Research Center for Nuclear Physics, Osaka University. 
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